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ABSTRACT 



We address the general problem of the luminosity-specific planetary nebula (PN) num- 
ber, better known as the "a" ratio, given by a = Np-^/L gll \, and its relationship with age 
and metallicity of the parent stellar population. Our analysis relies on population synthesis 
models, that account for simple stellar populations (SSPs), and more elaborated galaxy mod- 
els covering the full star-formation range of the different Hubble morphological types. This 
^■ig . theoretical framework is compared with the updated census of the PN population in Local 

' Group galaxies and external ellipticals in the Leo group, and the Virgo and Fornax clusters. 

\Q \ The main conclusions of our study can be summarized as follows: 

i) according to the Post-AGB stellar core mass, PN lifetime in a SSP is constrained by 
three relevant regimes, driven by the nuclear (M corc > 0.57 M©), dynamical (0.57 M© > 
M C ore > 0.55 Mo) and transition (0.55 M© > M coro > 0.52 M Q ) timescales. The lower 
limit for M coro also sets the minimum mass for stars to reach the AGB thermal-pulsing phase 
and experience the PN event; 

ii) mass loss is the crucial mechanism to constrain the value of a, through the definition of 
the initial-to-final mass relation (IFMR). The Reimers mass-loss parameterization, calibrated 
on Pop II stars of Galactic globular clusters, poorly reproduces the observed value of a in 
late-type galaxies, while a better fit is obtained using the empirical IFMR derived from white- 
dwarf observations in the Galaxy open clusters; 

rS \ Hi) the inferred PN lifetime for Local Group spirals and irregulars exceeds 10000 yr, 

■ which suggests that A/ corc < 0.65 M Q cores dominate, throughout; 

iv) the relative PN deficiency in elliptical galaxies, and the observed trend of a with 
galaxy optical colors support the presence of a prevailing fraction of low-mass cores (M core < 
0.55 Mq) in the PN distribution, and a reduced visibility timescale for the nebulae as a conse- 
quence of the increased AGB transition time. The stellar component with M core < 0.52 M©, 
which overrides the PN phase, could provide an enhanced contribution to hotter HB and Post- 
HB evolution, as directly observed in M 32 and the bulge of M 3 1 . This implies that the most 
UV-enhanced ellipticals should also display the lowest values of a, as confirmed by the Virgo 
cluster early-type galaxy population; 

v) any blue-straggler population, invoked as progenitor of the Af corc > 0.7 Mq PNe in 
order to preserve the constancy of the bright luminosity-function cut-off magnitude in ellip- 
ticals, must be confined to a small fraction (few percents at most) of the whole galaxy PN 
population. 

Key words: galaxies: evolution - galaxies: stellar content - galaxies: spiral - Galaxy: funda- 
mental parameters - ISM: lines and bands 



1 INTRODUCTION 

Diffuse intracluster luminosity lUson. Boughn. & Kuhnll99lh and 
other faint surface-brightness features detected at large distances 
from the center of isolated and cluster galaxies <Hui et allfl993l: 



Mihos et alj|2005l) may provide a valuable record of the mecha- 
nisms that led to the assembly and formation of the cosmic struc- 
tures at the different hierarchical scales. 

Stellar streams, like the case of Sgr and CMa in the Milky 
Way (e.g. iMartin et alJl2004 . and "runaway" halo stars iBlaauwl 
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ll96lUKeenan & Duftonlll983l : lAllen & KinmarjEool are, in this 
sense, excellent examples of the effect of galaxy tidal interac- 
tions and dynamical relaxation processes, that spread stars well 
outside the bright body of their parent systems. Furthermore, the 
evidence of quiescent on-going star formation in low-density en- 
vir onments, such as dwarf or very low surface-brightness galax- 
ies iHunter & GallaeheJl 9851: ICaldwellll 9951: Icellone & Buzzonil 
2005), indicates that the transition between coherent stellar systems 
and diffuse background might be not so sharp. 

In this framework, the study of Planetary Nebulae (PNe) is of 
special interest as these objects are efficient tracers of their under- 
lying parent stellar population even in those regions where the stel- 
lar plo t is too faint to be detected against the night-sky brightness 
I Arna boldi et al 12 003 \ Relying on standard narrow-band imagery, 
PNe ha ve been co nfidently detected out to 70 kpc from the center of 
Cen A iHui et al] 1 19931) and other bright galaxies in the Vir go and 
Fornax clusters, about 15 Mpc away l Arnaboldi et al. 2003, 2004; 
iFeldmeier et all2003lEo04l : IPeng et all2004l:lAguerri et all2005T) . 
On the other hand, optimized multi-slit imaging recently pushed 
this distance limit even farther, re aching the Coma cluster PNe, at 
about 100 Mpc distance jGerhard et all2005l) . 

One open question in this approach is the link between the size 
of the PN component and the sampled (bolometric) luminosity of 
its parent stellar population. This ratio, often referred to in the lit- 
erature as the "luminosity-specific PN number density" or, shortly, 
the "a ratio" ljacobvl ll980). provides the amount of light associ- 
ated to any observed PN sample, and it is the first assumption when 
using PNe as tracers of the spatial distribution and motions of the 
parent stars. Observationally, there is strong evidence for a to cor- 
relate with galaxy color, the reddest ellipticals being a f actor 5 to 
7 poorer in PNe per u nit galaxy luminosity than spirals ( Peimbert 
ll99rtlHuiet"aT]ll993l) . This trend is at odds with a nearly "uni- 
versal" PN luminosity function (PNLF), as observed for galaxies 
along the whole Hubble morphological sequence iCiardu llo et alJ 
l2002al) . 

Population synthesis models offer a useful reference tool, in 
this regard, as they allow us to probe the theoretical relationship of 
a with the different distinctive parameters of a stellar population, 
thus overcoming most of the observational uncertainties related to 
systematic bias, incomplete sampling etc.. In addition, a theoretical 
assessment of the a ratio does not require, in principle, any exact 
knowledge of the PN LF, still largely uncertain at its faint magni- 
tudes Ijacobv & De MarcofedbllljacobvfcOOd) . 

In this paper we want to study the luminosity-specific PN 
number density from a more general approach based on a new set 
of stellar population models. Our calculations rely on the Buzzoni 
ll989l hereafter B89) and lBuzzonil < 19951) synthesis code, includ- 
ing the recent extensi on to template galaxy models of composite 
star formation (Buzzoni 2002, 2005). Our predictions for the a ra- 
tio will also be compared with the observed PN population of Lo- 
cal Group galaxies, in order to validate our adopted scenario for 
Post-asymptotic giant branch (PAGB) evolution. We will also ad- 
dress the puzzling problem of the poorer PN population in elliptical 
galaxies iHui et all 19931) : this very interesting feature might be the 
signature of a phase transition related to horizontal branch (HB) 
evolution, with important consequences on the expected ultraviolet 
evolution of galaxies and other stellar systems in metal-rich envi- 
ronments. 

Our paper will be organized as follows: in Sec. |2|we shall 
present the simple stellar population (SSP) theory, and discuss the 
relevant parameters that constrain a. Section [3] is devoted to the 
expected theoretical values of the a parameter for galaxies of dif- 



ferent morphological type. In Sec. |4| we compare the theoretical 
results with observations of the PN population in Local Group and 
external galaxies in the Leo group, the Virgo and Fornax clusters. 
In Sec.|5| we study the correlation between a and other spectropho- 
tometric indices for early-type galaxy diagnostic, like the Lick Mg2 
index and the (1550 — V) ultraviolet color. Finally, in Sec.|6]we 
summarize the constraints obtained for the PAGB evolution and the 
relevant timescale for PN evolution. 



2 SIMPLE STELLAR POPULATION THEORY 

The SSP theory developed by Re nzini & Buzzonil ll986l hereafter 
RB86) and B89 allows us to compute the value of a for a coeval 
and chemically homogeneous stellar aggregate. Assuming a stan- 
dard power-law IMF such as N* (M,) oc M^ s dM*, the expected 
number of stars that populate the "j"-th Post-main sequence (PMS) 
phase, of duration Tj, can be written as: 



Nj = AM T q | A/ TO | Tj 



(1) 



where Mto is the time derivative of the MS turn off (TO) mass. 
The scaling factor A in eq. Q accounts for the total mass of the 
SSP (cf. B89), and relates to the total luminosity of the stellar pop- 
ulation, so that Nj oc Ltot Tj . 

When this approach is applied to the PN phase, it becomes: 



A?PN = B L to t TpN, 



(2) 



where B = AM, 



. .|M xo| /I/tot is the so called "specific evolu- 
tionary flux" (see RB86 and B89), and r PN is the PN visibility life- 
time, i.e. the time for the nebula to be detectable in [Om] and/or 
Ha surveys. The value of the "luminosity-specific PN number" a 
simply derives from eq. J2j as 

iVpN 



= Sr PN . 



(3) 



Therefore two parameters, that is B and t P n, set the value of the 
luminosity-specific PN number for a given SSP. 



2.1 The "specific evolutionary flux" B 

This parameter links the PN rate to the evolutionary properties 
of the parent stellar population. It contains, in fact, i) the phys- 
ical "clock" of the population (as both Mto and Mto depend 
on time), and ii ) the IMF dependence, because the ratio A/Ltot 
closely scales with the SSP M/L ratio, thus responding to the mass 
distribution of the comp osing stars. 

Following [5uzzonj| ll998l) . we can further arrange our defini- 
tion of B and write 



B 



AM-^\Mto\ 



£ 



(4) 



Ltot PMS fuel 

The r.h. side of eq. J4) derives from the RB86 "Fuel consumption 
theorem" and identifies the PMS contribution to the total SSP lu- 
minosity (C = I/pMs/^tot), as well as the absolute nuclear fuel 
spent during PMS evolution by stars with M» = Mto- 

As pointed out by RB86, B does not change much with metal- 
licity or IMF slope; this is shown in Fig.0 where we display the 
specific evolutionary flux, according to B89 SSP models, for dif- 
ferent IMF slopes and metallicity. One can see that models with 
Z = 1/20 to 2 Z Q and IMF slopes from 1.35 to 3.35 changes B by 
only a factor of two, at most, over a wide range of SSP age, from 1 
to 15 Gyr. 
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Figure 1. Specific evolutionary flux B, from eq. |4j, for B89 SSP models. 
Different metallicity sets (between Z ~ 1/20 and 2 Zq) are overplotted. In 
addition to the Salpeter case (s = 2.35), other IMF power-law coefficients 
are explored, as labeled on the plot. The value of B is given in units of 



The dependence of a on the IMF slope is through the change 
in £, i.e. the relative PMS luminosity contribution in eq. In 
general, a giant-dominated SSP will display a larger value of a be- 
cause giant stars become more important in a "flatter" IMF and £ 
tends to increase, compared to the Salpeter s — 2.35 case. On the 
other hand, a younger age acts in the sense of decreasing the PN 
rate per unit SSP luminosity, because the MS luminosity contribu- 
tion increases, the absolute amount of "fuel" burnt by PMS stars 
increases (cf., e.g. Fig. 3 in RB86), and £ decreases. 



2.2 The PN lifetime t pn 

This quantity depends both on the chemio-dynamical properties of 
the ejected material, giving rise to the nebular envelope, and on the 
stellar core-mass evolution, which is responsible for the "firing up" 
of the nebular gas. 

The first self-con sistent the o retical model of PN core 
evolution is due toJPacz^nskil ^^1), then f ollowed by 
contrib uti ons from |§ chonbernej ^98lL |l983|). [ G orn v et alJ 
Jl994l). IVassiliadis & Wool T\9<m . IStanghellinil <1995l) 
Stanahellinj^^^^^^ (2000), and more recently bv lMarigo^Htl] 
(2001) and IPerinotto et alJ <20043) . Theory predicts that the PN 
properties strongly depend on the core-mass distribution of PAGB 
remnants, the latter being the result of the initial-to-final mass 
relation (hereafter IFMR), as modulated by metallicity and mass 
loss. The so-called "transition time" for PAGB stars to reach 
the high-temperature region (T e fi > 30 000 K), after the nebula 
ejection, may also play a role on the emission properties and the 
timescale of the PN event. 

A simple estimate of PN core lifetime can be derived from 
the energetic budget available to PAGB stars. This is co mputed 
in Fig . 13 based o n a m odel collection fr om Paczvhski 1 1971), 
SchonberneJ 1 1 98 ill 1 983 ) andlVassiliadis & WooJ<1994 Indeed, 
one sees that lPaczvriski'; ll97ll) original results consistently meet 
the more recent and sop histicated stellar tracks, like those of 
IVassiliadis & Woocl 1 19941) . that account for stellar envelope ejec- 
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Figure 2. Theoretical fuel consumpt ion for stars alon g the PAGB evolu- 
tion according to different model sets: Paczvriski 1 1971, pentagon markers), 
ISchonb ernei 1 1981, 1983, squares and rhombs)|Vassiliadis & Wood 1 19941 
dots and triangles). The different markers for the same model source refer 
to the prevailing case of a H or He thermal pulse terminating the AGB evo- 
lution, as labeled. Fuel is expressed in Hydrogen-equivalent solar mass, i.e. 
1 g of H-equivalent mass = 6 10 18 ergs (cf. RB86), and a solar metallicity 
is assumed in the models. A smooth analytical relation matching the data, 
according to eq. |5j is plotted as a solid curve. 



tion either in the case of a He or H thermal pulse at the tip of AGB 
evolution. 

From B89, an analytical function that reproduces the PAGB 
fuel consumed as a function of the core mass for the different the- 
oretical data sets of Fig.|2|is 

-5.22 



PAGB fuel = (M core /0.163)" 



\H M, 



oj 



(5) 



where M CO rc is in solar unit and fuel in Hydrogen-equivalent so- 
lar mass (cf. Fig. [2] for details). The PAGB stellar lifetime can be 
defined as 



tpagb ^ (PAGB fuel/Ws) 



(6) 



being ^pagb the luminos ity of the stellar core at the onset of the 
nebula ejection (see e.g. |Paczvnskilll97ll) . In general, tpagb is 
shorter than the estimated dynamical timescale for the nebula evap- 
oration, which is about Td yn ~ 30 000 yr (cf. SchonberneJ ll983t 
|Phillipslll989l) . and this implies that for SSPs of young and inter- 
mediate age r PN ^ tpagb- 

On the other hand, when the SSP age increases, this assump- 
tion may no longer hold, and at some point, when PAGB stellar 
core mass decreases and M corc < 0.57 Mq, tpagb starts to ex- 
ceed Td yn . Henceforth, the PN evolution is driven by the dynamical 
timescale, reduced by the transition time. 



2.3 PAGB core mass and PN evolution 

As M C orc constrains both the fuel and PAGB luminosity, from 
which tpagb derives, 1 it is clear that mass loss mechanisms, at 
work along the giant-branch evolution, set the leading parameters 



1 In force of eq. |5}, and also considering that <?pagb oc A7 C orc in eq. 
lPaczvriskll97lh . we have that tp^qb oc Mcotc 



4 A. Buzzoni, M. Arnaboldi & R.L.M. Corradi 




0.5 



10 
[Gyr] 



5 10 
Age [Gyr] 



5 10 
Age [Gyr] 



Figure 3. Time evolution of the stellar mass at 
some tipping points across the H-R diagram for 
SSPs of different metallicity Z, about the so- 
lar value, as labeled in each panel. Upper strips 
in each panel are the theoretical loci for stellar 
mass at the tip of the RGB evolution, ( Mrgb ) 
according to a Reimers mass loss parameter in 
the range 0.3 ^ T] ^ 0.5, as labeled in the mid- 
dle panel, for general reference. Lower strips 
mark the locus for stellar mass at the onset of 
PAGB evolution (AfpAGB), again for the same 
reported range of the mass loss parameter r\. 
The minimum mass for stars to reach the AGB 
thermal pulsing phase (and eventually produce 
a PN) is marked by the dashed line, accord ing 
to lDorman et"ai] <1993l) and lBlockeJ ll995t) . 



for PN evolution. In this sense, it is of paramount importance to es- 
tablish a suitable relationship between the initial (i.e. Mi = Mto) 
and final (Mf = Mpagb) mass of stars along the whole SSP evo- 
lution. 

A firm settlement of the IFMR is a long-standing problem, that 
has been addressed both theoretically (Iben & Renzini 1983, here- 
after IR83, Dominguez et al. 1999, Girardi et al. 2000) and obser- 
vationally iweidemann & Koestet 1 1 983t IWeidemanrjl 19871 12000L 
hereafter W00; IClaver et all200ltlKalirai et all2005t) . 



2.3.1 Mass loss and M corc : the theoretical approach 

A dire ct evaluation of the mass loss, according to the iReimerJ 
theoretical parameterization, is shown in Fig. [5] In this fig- 
ure, we plot the expected value of stellar masses at some tipping 
points of SSP evolution, for three relevant values of metallicity 
around the solar value. In particular, A4rgb is the mass of stars at 
the end of red giant branch (RGB) evolution, which occurs after the 
first important mass-loss episode experienced by stars at the low- 
gravity low-temperature regime; Mpagb is then the mass of the 
stars leaving the AGB, after the second stronger mass-loss episode. 
The trend of both quantities is tracked vs. SSP age, for a range of 
the Reimers mass-loss parameter, r\, between 0.3 and 0.5. 2 

Along the SSP evolution shown in Fig.|5| a substantial fraction 
(up to 50%) of the stellar mass is lost during the AGB phase (cf. the 
difference AM, — A/rgb — Mpagb on the plots); furthermore, 
with increasing age, the core mass of PAGB stars approaches (or 
even crosses) the limit of Mpagb — 0.52 Mq, which is the mini- 
mum mass required by models for stars to experience the so-calle d 
AGB "thermal-pulsing" phase JPorman et alJl993tlBiockeJl995l) . 
During this phase, stars venture in the region of Mira variables 
and end their AGB evolution with a quick envelope ejection, likely 
driven by dynamical instability (the so-called "superwind phase" of 
lRenziniH 98ll) . and the subsequent PN stage (see IPaczvnskil 1 97oL 
and I R831 for an exhaustive discussion of the process and its vari- 
ants). We shall discuss the implication of an inhibited AGB phase 
on PN evolution in what follows. 



2 A value of r\ ~ 0.4 ± 0.1 is typically required in order to reproduce the 
observed colour-magnitude (c-m) diagram of old Galactic globular clusters 



2.3.2 Mass loss and M corc : the empirical approach 

The IFMR can be derived empirically from the mass estimate of ob- 
served white dwarfs in nearby open clusters (like Hyades or Prae- 
sepe), coupled with the value of Mto and cluster age, as obtained 
from isochrone fit of the cluster c-m diagram. 

Quoting WOOL there is evidence that observed white dwarf 
masses, for low- and intermediate-mass stars, "coincide almost ex- 
actly with the new theoretical predictions of the core masses at 
the beginning of the thermal pulsing AGB" (A/tp). In addition, 
"it is presumed and supported weakly by the empirical data that 
this closeness of the final mass to the first-thermal pulse core mass 
relation continues also to higher initial masses", possibly up to the 
limit of SN onset (roughly placed about 7 Mq). 

The first clai m of t his analysis is confirmed by Fig.|4] where 
we compare the WOO IFMR and the theoretical estimates of 
Mtp from an updated set of stellar tracks for Pop I stars by 
Waaenhuber & Weiss 1 1994) and from the original analytical re- 
lation by IR83 for intermediate-mass stars, namely 



M T p = 0.59 + 0.0526 A'/;. 



(7) 



It is also clear from Fig. [4] that, for a standard range of 
Re imers mass-loss parameters fitted to Galactic globular clusters, 
the IR83 theoretical IFMR overestimates the value of Mf for young 
(t < 2 Gyr) SSPs, requiring a value of 77 > 1 to match the lWOfJ 
empirical relation. 

Given the discrepancy between the predicted M core based on 
the standa rd ma ss loss theory a /a lReimersI 1 19751) and the obser- 
vations by WOO, we compute the t pagb values according to a) a 
theoretical IFMR according to B89 SSP models with r\ = 0.3, ex- 
tended to younger ages through the IR83 relation 



M c . 



0.53 7/ 



-0.082 . n , r -0.35/ a , 

+ 0.15?? (Mto 



(8) 



and h) the WOO empirical IFMR, as displayed in Fig. [4] 



iFusi Pecci & Renzini 197(1) . 



2.3.3 A critical core mass range: 0.52 Mq ^ M core ^ 0.55 Mq 

When the core nuclear lifetime ceases to be the driving parame- 
ter for PN visibility (that is for M corc < 0.57 Mq), we have 
that tpn — (Vdyn — Ttt), where Ttt is the transition time of the 
stellar core to become hot enough such as to "fire up" the nebula. 
Models show that, for M corc < 0.55 Mq, Ttt abruptly increases 
from its typical value range of 200-2000 yr up to exceeding Td yn 
lSch6nberneJl983l : IVassiliadis & Wooo1l994l) . In this case, by the 
time the stellar core is there ready to heat, the gaseous shell has al- 
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Figure 4. The initial-to-final mass relation according to different calibra- 
tions. The solid strip is the theoretical relation of IR83 for a standard mass 
loss parameter r\ in the range between 0.3 and 0.5, as labeled on the plot. 
Small dots report the individual values as from the |b89 SSP models of 
Table HI and the same Reimers parameters. Short- and long-dashed curves 
are the theoretical loci for star s to set on the AGB thermal pulsing phase 
(M TP ), according to llR83l and lwagenhuber & WeisJ fl994l> (WW94V Fi- 
nally, big dots and solid curve report the WOO empirical relation based on 
the mass estimate of white dwarfs in Galactic open cluster. 



ready evaporated, thus preventing the nebula ignition. The dynam- 
ical timescale itself is slightly influenced by mass loss as a slower 
envelope expansion, driven by radiation pressure, is expected if r\ 
incr eases such as to terminate AGB evoluti on at l ower luminos- 
ity JVassiliadis & Woodlll994l iMarigo et aljEoOll Ivillaver et"ai] 
120021) . 

In general, this scenario leads to conclude that in old SSPs 
PN visibility might be greatly reduced or even fully inhibited as a 
consequence of a delayed ho t-PAGB evolution of the stellar core 
IStanghellini & Renzinl2000l) . 

In addition, one has also to account for a further critical thresh- 
old in the evolutionary framework when stars escape the AGB 
thermal-pulsing phase, for M core < 0.52 Mq. The lack of a 
full AGB development leads to a ra nge of Post-HB ev olution- 
ary paths, 3 as discussed in detail by iGreggio & Renzirill il990h . 
One relevant case, in this regard, is that of "AGB-manque" stars, 
that directly set on the high-temperature white-dwarf cooling se- 
quence after leaving the HB, thus missing, partially or in toto, 
the AGB phase. First important hints of this non-standard evolu- 
tionary framework can be found in the original work of Gingold 
19743, and a number of later contributions tried to better as- 
sess the critical parameters (m ass los s in primis) involved in the 
process iCastellani et alJll992L Il995t l2006tjDorman et alJll993t 
iD'Cruz et alll996tlYi. Demarque. & OemleJl998h . 

The impact of this composite scenario on the PN forma- 
tion mechanisms is probably still to be fully understood; it seems 
clear, however, that a full completion of AGB evolution, up to the 
thermal-pulsing luminosity range, i s the most vi able step to culmi- 
nate with the nebula event (see also Kwok 1994, for a more elabo- 
rated discussion in this regard). 



3 From the physical point of view, this would correspond to the He+H 
double-shell burning regime for low- and intermediate-mass stars. 
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Figure 5. The luminosity-specific PN number for SSP models of Table HI 
(both for r) = 0.3 and 0.5 and different metallicity, as reported top left) 
compared to the PAGB stellar core mass. Overplotted are also the expected 
calibration assuming the theoretical IFMR of IR83 and the empirical one 
from [WOO, as labeled. Note the clean relationship in place, with M COIe 
being the leading parameter to constrain a. 
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Figure 6. Theoretical time evolution of the luminosity-specific PN number 
for SSP models of Table HI tboth for rj = 0.3 and 0.5 and the different 
metallicity values, as labeled top left on the pl ot) com pared to the expected 
calibrations assuming the theoretical IFMR of llR83l and the empirical one 
from lWOOl 

2.4 The luminosity-specific PN number in Simple Stellar 
Populations 

According to the previous discussion, the luminosity-specific PN 
number in a SSP can eventually be written as 

_ J Srnin[r P AGB, (r dyn - r u )], if Ttt < T dyn 

y otherwise. 

In general, r c i yn sets a safe upper limit to the a value. For a 
Salpeter IMF, and from the data of Fig.pJ we derive 

a max ~ 1.8 nr 11 x 30 000 = - J™ . (io) 

1.85 10 b Lq 

This value is virtually independent from metallicity. 
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Table 1. Luminosity-specific PN number for Salpeter SSPs( a ' 
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-8.16 


-7.22 


1.0 










-7.98 










-7.60 


-6.86 


2.0 




-7.29 


-7.41 


-7.41 


-7.45 




-6.99 


-7.02 


-7.08 


-7.11 


-6.60 


3.0 




-7.08 


-7.15 


-7.18 


-7.17 




-6.70 


-6.75 


-6.76 


-6.85 


-6.47 


4.0 




-6.90 


-6.94 


-6.94 


-6.98 




-6.53 


-6.56 


-6.55 


-6.68 


-6.38 


5.0 




-6.74 


-6.85 


-6.88 


-6.84 




-6.39 


-6.40 


-6.42 


-6.54 


-6.32 


6.0 




-6.64 


-6.72 


-6.73 


-6.73 




-6.29 


-6.30 


-6.30 


-6.45 


-6.30 


8.0 


-6.47 


-6.46 


-6.52 


-6.57 


-6.56 


-6.29 


^-6.28 


sC-6.28 


^-6.29 


-6.28 


-6.27 


10.0 


-6.34 


-6.31 


-6.39 


-6.45 


-6.44 


^-6.27 


1 






-6.27 


-6.27 


12.5 


-6.29 


-6.28 


-6.29 


-6.30 


-6.31 






no PNeW 




-6.27 


-6.27 


15.0 


^-6.28 


^-6.27 


^-6.28 


^-6.29 


-6.27 


1 








^-6.27 


^-6.27 



( a ) The listed quantity is loga = logJVpN — logC-LsSp/^©)' 

W Planetary nebulae not expected to form for these age/metallicity combinations as AfpAGB Ss 0.52 Mq. 



The theoretical luminosity-specific PN number for B89 SSPs 
of different age, metallicity and mass-loss parameter is summarized 
in Fig. [5] and [6] by matching the IFMR prescriptions according 
both to case (a) and (b), as in previous discussion. As expected, 
Fig.0shows that M CO rc is indeed the leading parameter constr ain- 
ing a ; note however the relevant difference between the IR83 and 
WOO models, with the latter reaching a fixed core mass at younger 
age, w hich corresponds to a brighter SSP luminosity. Compared to 
IR83, therefore, the WOO IFMR predicts in general a lower value 
of a for fixed value o/M coro . 

In this framework, the role of other SSP distinctive parame- 
ters, like metallicity and mass loss, only enters at a later crucial 
stage of the analysis, when tuning up the clock that links M corc to 
SSP age. This is displayed in Fig.|6| where the theoretical evolution 
of a is traced for SSPs spanning the full range of metal abundances 
and Rei mers r \ parameter. Again, fhe|TR83 r\ = 0.3 and 0.5 cases 
and the IWOOl IFMR are compared in the figure. Whe n the varia- 
tion of a is predicted as a function of the SSP age, the WOO model 
leads to a systematically higher PN number per unit SSP luminos- 
ity, compared to the IR83 case, due to a lower value of PAGB core 
mass assumed. 

A summary of our calculations is reported in Table Q A 
Salpeter IMF was adopted in the models, but a change in the power- 
law index with respect to the canonical value of s = 2.35 can easily 
be accounted for, following eq. J4j: 



A log a — log 



a B 

«Sal 



(ID 



This implies that, for a giant-dominated SSP (s = 1.35), loga is 
~ 0.04 dex higher, and the opposite happens for a dwarf-dominated 
SSP (s = 3.35), for which log a is ~ 0.10 dex lower than for the 
Salpeter case (see Fig.0. 

The effect of enhanced mass loss on the Fig. |6| models can 
be quantified in roughly A log a ~ +0.4 dex for a change of 77 
from 0.3 to 0.5 and fixed SSP age. This is a consequence of a lower 
core mass and a correspondingly slower PAGB evolution (see foot- 
note0, recalling that A log a — A log tpagb- 

A different value for the evaporation timescale directly reflects 
on amax. From eq. 1101 . as A log a max = A log Td yn , one has for 
instance that Aloga max ~ —0.18 dex when Ta yn is decreased, 
say, to 20000yr. 



Given the extreme uncertainty of theory in quantitatively as- 
sessing the PAGB transition time, in our SSP models we do not 
explicitly account for the quick drop of luminosity-specific PN 



< 0.55 M s PAGB core 



number for the 0.52 M Q < M c , 
mass range. Operationally, as far as the PN evolutionary regime 
is driven by Td yn , we neglect the r t t contribution and simply as- 
sume q = Qfmax, 4 while for Afcorc ^ 0.52 Mq no PN events are 
expected to occur at all and a = 0. This crude simplification has 
negligible consequences when modeling late-type galaxies, while 
the evolution in the 0.52 < M core < 0.55 Mq mass range be- 
comes important for early-type galaxy models, as we shall discuss 
in more detail in Sec. [5] 

As a concluding remark, we must also recall that a vanishes 
for SSP ages t < 10 8 yr, when the prevailing high-mass stars of 
5-7 Mq or higher override the PN phase and end up their evolution 
as Supernovae. 



3 THE LUMINOSITY-SPECIFIC PN NUMBER IN 
GALAXIES 

The evolutionary properties of PNe in SSPs are the basis for ex- 
tending the analysis to a wider range of star formation histories, 
as it happens in real galaxies . We u se the template galaxy models 
developed by Buzzoni 1 2002, 2005), which ensure a self-consistent 
treatment of the PN evolution and the photometric properties of the 
stellar populations in the parent galaxy. In this framework, colors 
and morphological features along the Hubble sequence for early- 
and late-type systems were reproduced tracing the individual lu- 
minosity contribution for the bulge, disk and halo component (see 
lBuzzorj200ll2005l for additional details) . 

Following Sec. 12.41 the PN evolution has been implemented 
in the models in a semi-analytical way, by fitting the SSP data in 
TableQJfor the r) = 0.3 and 0.5 cases, and for the WOO empirical 
IFMR. The whole data grid can be fitted within a 6% formal un- 
certainty in the value of a [i.e. er(loga) = ±0.025 dex] over the 
age/metallicity range by 



In this case, an upper limit is marked for a in TableQ 
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Table 2. The luminosity-specific PN number for template galaxy models'™' 



Age 


M bo l 


B-V 


log a 3 


log a w 


M b ol 


B-V 


log a 3 


\oga w 


Mbol 


B-V 


log a 3 




[Gyr] 
































F. 








Sn 








Sh 




1.0 


-23.10 


0.66 


-7.84 


-6.87 


-23.22 


0.58 


-7.94 


-6.98 


-23.11 


0.55 


-7.99 


-7.02 


2.0 


-22.48 


0.72 


-7.40 


-6.56 


-22.67 


0.63 


-7.51 


-6.67 


-22.66 


0.58 


-7.59 


-6.74 


3.0 


-22.13 


0.76 


-7.16 


-6.43 


-22.36 


0.65 


-7.27 


-6.54 


-22.42 


0.59 


-7.36 


-6.63 


4.0 


-21.89 


0.79 


-6.98 


-6.36 


-22.14 


0.67 


-7.09 


-6.47 


-22.25 


0.60 


-7.20 


-6.57 


5.0 


-21.70 


0.81 


-6.84 


-6.32 


-21.97 


0.69 


-6.95 


-6.43 


-22.13 


0.61 


-7.08 


-6.54 


6.0 


-21.54 


0.83 


-6.73 


-6.30 


-21.84 


0.70 


-6.84 


-6.41 


-22.03 


0.62 


-6.97 


-6.52 


8.0 


-21.30 


0.86 


-6.55 


-6.28 


-21.62 


0.72 


-6.66 


-6.38 


-21.88 


0.63 


-6.81 


-6.49 


10.0 


-21.11 


0.88 


-6.41 


-6.27 


-21.45 


0.74 


-6.53 


-6.37 


-21.76 


0.63 


-6.68 


-6.48 


12.5 


-20.92 


0.90 


-6.29 


-6.28 


-21.29 


0.75 


-6.42 


-6.37 


-21.65 


0.64 


-6.58 


-6.48 


15.0 


-20.76 


0.92 


-6.29 


-6.29 


-21.15 


0.76 


-6.41 


-6.38 


-21.57 


0.65 


-6.57 


-6.49 
















Sri 








Tm 




1.0 


-22.80 


0.55 


-7.99 


-7.03 


-22.12 


0.51 


-8.06 


-7.09 


-20.29 


0.30 


-9.03 


-7.73 


2.0 


-22.47 


0.55 


-7.65 


-6.80 


-22.05 


0.48 


-7.80 


-6.94 


-21.06 


0.34 


-8.57 


-7.55 


3.0 


-22.34 


0.54 


-7.46 


-6.72 


-22.11 


0.46 


-7.66 


-6.91 


-21.49 


0.36 


-8.31 


-7.41 


4.0 


-22.26 


0.54 


-7.33 


-6.68 


-22.18 


0.46 


-7.56 


-6.89 


-21.80 


0.38 


-8.13 


-7.30 


5.0 


-22.22 


0.54 


-7.22 


-6.66 


-22.25 


0.46 


-7.48 


-6.88 


-22.04 


0.39 


-7.98 


-7.21 


6.0 


-22.19 


0.54 


-7.14 


-6.65 


-22.32 


0.46 


-7.41 


-6.87 


-22.24 


0.40 


-7.86 


-7.15 


8.0 


-22.15 


0.54 


-7.00 


-6.64 


-22.44 


0.47 


-7.29 


-6.84 


-22.55 


0.42 


-7.67 


-7.05 


10.0 


-22.13 


0.55 


-6.90 


-6.63 


-22.54 


0.48 


-7.19 


-6.82 


-22.79 


0.43 


-7.53 


-6.99 


12.5 


-22.12 


0.55 


-6.81 


-6.63 


-22.64 


0.49 


-7.10 


-6.79 


-23.03 


0.45 


-7.39 


-6.93 


15.0 


-22.11 


0.56 


-6.79 


-6.63 


-22.73 


0.49 


-7.05 


-6.77 


-23.23 


0.46 


-7.29 


-6.88 



( a ' Models are for a Salpeter IMF; 

log c*o3 = luminosity-specific PN number assuming a theoretical IFMR with a fixed Reimer s mas s loss p arameter r) = 0.3; 
log otyy = luminosity-specific PN number assuming an empirical IFMR according to Weidemann 1 2000). 



log a' = (1.52 - 0.07 z) log tg - 0.07 log 2 (3.4 z) 
-0.1/tg + 2.0 (77- 0.3) -7.80, 



(12) 



where tg is the SSP age in Gyr, an d z = Z/Zq . For solar metallic- 
ity, eq. 1121 reproduces the IR83 calibration over the extended age 
range of TableQ Similarly, a fit to the WOO relation is: 



log a = -0.6 (log tg - \f - 6.27. 
In addition, we also assume that 



min [a , a„ 




for Tig ^0.1 
iortg < 0.1 



(13) 



(14) 



where a max is given by eq. 1 1 ()t . 

The absolute number of PNe in a SSP of total (bolometric) 
luminosity Lssp simply becomes A^pn (SSP) = a Lssp , and for a 
star-forming galaxy of age t we write 



M»n(*)= / a(r)L SS p(r) SFR{t-r)dr. 
Jo 



(15) 



Therefore, the global luminosity-specific PN number for the 
galaxy stellar population can be computed as 



a{t) 



■A/pn(«) _ So a ( T ) issp(r) SFR{t - 



dr 



Jgal 



J Lssp (t) SFR(t-r) dr 



(16) 



The results for the whole Hubble sequence, from type E to Im, 
are summarized in Fig.|7jand Table|2|for the 77 = 0.3 case (003 in 
the table) and the IWOOl IFMR {aw)- In Tabled we provide also 
the integrated B — V color for the parent galaxy, and the absolute 
bolometric luminosity of the system, assuming a total stellar mass 



of Mgai = 10 U M© at 15 Gyr (see lBuzzonill2005L for a detailed 
definition of this quantity). 5 

The two evolutionary scenarios produce large differences in 
a(t) for the Hubble morphological types: the WOO model pre- 
dicts larger PN populations and a shallower time dependence for 
a than for the 77 = 0.3 case. The effects are larger for late-type 
galaxy models, where the impact of the different IFMR details on 
intermediate- and high-mass stars is stronger. Also, the nucleated 
late-type galaxies (types Sa => Sd) approach the evolutio n of el- 
liptica ls at early epochs. As discussed in some detail in iBuzzonil 
(2005), this effect is due to the prevailing photometric contribution 
of the bulge stellar component, when t—*0. 

The general prediction from these models is that a is expected 
to decrease in young and/or star-forming galaxies, compared to 
more "quiescent" early-type systems as a consequence of a smaller 
population of PNe embedded in a higher galaxy luminosity per unit 
mass (i.e. a lower M/L ratio). 



5 The lBuzzoni <2005l) galaxy models are computed with a standard mass- 
loss parameter r) = 0.3, and the integrated magnitudes and colors reported 
in Table |2) refer to this case. The effect on galaxy B — V by adopting the 
alternative WOO IFMR can be estimated in A(B— V) ~ -0.02 mag, that is 
w ith a little shift toward bluer colors (see e.g., the calibration from Fig. 28 
in lBuzzoni[l995T) . The value reported in Table l2l for loga^- should also 
be increased by an additional A log«vf — +0.04 dex if one takes into 
account slightly fainter galaxies (some -10% in bolometric luminosity) in 
the WOO framework. 
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1 5 10 15 i 5 10 15 

Age [Gyr] Age [Gyr ] 

Figure 7. Theoretical time evolution of the luminosity specific PN density (a) for the Buzzoni 1 2005) template galaxy models along the whole E-Sa-Sb-Sc- 
Sd-Im Hubble morphological sequence. Models in left panel assume an IFMR as from the standard mass loss parameter r\ = 0.3, while those in the right 
panel rely on the empirical relation from WOO. Note, in the latter case, the much shallower evolution of a. In the two panels, bulge-dominated spirals tend to 
approach the evolution of ellipticals at early epochs due to the increasing bulge contribution to the global galaxy luminosity. 



3.1 Comparing with the observations: PN luminosity 
function and completeness corrections 

The empirical evidence, from the 5007 A [Olll] PN luminosity dis- 
tribution, indicates a nearly constant value for the bright cut-off 
magnitude (A/*) of the PNLF. This feature is actually found to 
be near l y inva riant with galaxy type and age and, as suggested by 
Ijacobvl ll989l) . can be effectively used as a standard candle to de- 
termine extragalactic distances. 

When the emission-line fluxes are converted into equivalent V 
magnitudes via the formula 



m [0 in] = -2.5 log F [0 in] - 13.74 



(17) 



ljacobvllT989l). the PNLF takes the shape of a double-exponential 
function ICiardulloet ail 1989) of the form 



log N(M) = 0.133 M + log [1 



e 3(M *" M) ]+ const, (18) 



with the bright cut-off magnitude placed at M*=— 4.47 mag plus a 
little metallicity correction, that scales with PN Oxygen abundance, 
such as 



AM* = 0.928 [O/H] 2 + 0.225 [O/H] + 0.014 



(19) 

lDopitaetalll992l) . After correction, the inferred PN distances in 
external galaxies are consistent within 0. 1 dex with those obtained 
using the Cephei ds method for a large o bserved sample of galaxies 
and galaxy types ICiardullo et alj2002^ . 

There is no theoretical explanation of this "universal" prop- 
erty of the PN distribution, although it has been questioned (e.g. 
Mend ezet all 19931) that it might depend on the sample size, being 
therefore a mere statistical effect. Furthermore, on the theoretical 
side, some change of M* with the age of the PN parent stellar pop- 
ulation should be expected, its precise amount strongly depending, 
howe ver, on the model assumptions (see especially Mar igo et alj 
120041 for the most striking results on this line). 

While the bright end of the PNLF has been studied in a 
large number of galaxies, not much is known about the PNLF 
shape at fainter magnitudes, as most Local Group (LG) surveys 
are complete down to 4 ma g from the bright cut-off luminos- 
ity ICorradi & M agrini 2006), and this limit becomes obviously 



1 

0.8 
0.6 
0.4 
0.2 



10 - 



100 



I 1 1 1 I 1 1 1 I 
Cumulative 

Fraction 




Figure 8. Upper panel: the cumulative fraction of PNe in the different mag- 
nitude bins with respect to the luminosity-function bright cut-off (M*) for 
the double-exponential fit of the PNLF, as in eq. j 1 81 (solid curve) and for 
the empirical SMC luminosity function according to Jacobv 1 2006) (dashed 
curve). Lower panel: completeness factor (CF = Wtot/^M— M*)) 
for the same c alibrat ions as in the upper panel. Also reported are the 
ICiardulloet all 1 19891) data for M 31 (star markers), and the relevant cor- 
rection factor for the 012.5 parameter. For better convenience, data are also 
listed in Tablel3l 



brighter for more distant systems. This implies a still large uncer- 
tainty in the extrapolation from the observed number of PNe to the 
whole PN population size, through eq. 1 1 81 . In particular, there is 
evidence of a dip in the PNLF of the SMC Ijacobv & De Marcel 
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Table 3. The PN Luminosity Functions 



M-M* 


Cumulative Fraction 


Completeness Factor 


[mag] 


Standard 


SMC 


Standard 


SMC 


u.u 








CO 


CO 


U.J 


0.009 


0.013 


113. 


78.3 


1 


0.025 


0.042 


39.8 


23.7 


| ^ 


0.046 


0.097 


21.6 


10.3 


2.0 


0.071 


0.155 


14.1 


6.46 


2.5 


0.100 


0.193 


9.99 


5.19 


3.0 


0.134 


0.225 


7.46 


4.44 


3.5 


0.173 


0.269 


5.76 


3.72 


4.0 


0.219 


0.310 


4.56 


3.23 


4.5 


0.273 


0.336 


3.66 


2.98 


5.0 


0.336 


0.366 


2.98 


2.73 


5.5 


0.409 


0.427 


2.45 


2.34 


6.0 


0.494 


0.535 


2.03 


1.87 


6.5 


0.593 


0.692 


1.69 


1.44 


7.0 


0.708 


0.845 


1.41 


1.18 


7.5 


0.843 


0.941 


1.19 


1.06 


8.0 


1 


1 


1.00 


1.00 



120021) and M33 <Magrini et all2000lJCiardullo et all20 04) at 4 and 
2.5 mag, respectively, below the bright cut-off (cf. Fig. [8] dashed 
curve in the bottom panel), but this is not observed in the M3 1 bulge 
ICiardullo etai] |2002a). The presence of the dip might depend 
on the galaxy star formation history (through the age distribution 
of PN progenitors), possibly witnessing the presence of a s ignifi- 
cant c omponent of young stars ICiardullo et all2004lMarigo et alJ 
120041) 

According to ICiardullo etaD <1989l) and Ijacofrd <198(J) . the 
faint-end tail of the standard PN LF agrees with the theoretical lu- 
minosity function o flHenize & Westerlundl iT963). in which a PN is 
modelled as a uniformly expanding homogeneous gas sphere ion- 
ized by a non-evolving central star. The number of nebulae in each 
luminosity interval should then be proportional to the PN lifetime 
sp ent within that luminosi t y bin. For a global PN lifetime of 30 000 
yr, iHenize & Westerlund 1 19631) predict that faintest PNe should 
locate about 8 magnitudes below M*. 

Deep observations in the SMC IJacobv & De MarcdEool 
|jacobvll2006l) . reaching more than 6 mag below the bright cut-off, 
show a significant decline of the number of PNe compared to what 
predicted by the double-exponential formula. Clearly, a more ex- 
tended sample from other nearby galaxies at comparable magnitude 
depth would be required to better assess this important problem. 
Unfortunately, in case of late-type galaxies, the [Olll] PN detec- 
tion is not quite a simple task; due to the ongoing star formation, 
only a very small fraction of the [Olll] emission-line sources in 
spirals and irregulars is represented by genuine PNe, in most cases 
counts are affected by Hll regions and supernova remnants. Identi- 
fication of PNe in LG galaxies requires additional constraints, such 
as point-like on-line emission combined with a non-detection in 
the off-line continuum, and a stronger [Olll] emission compared to 
Ha and/or [Nil] narrow-band luminosity. 6 In addition, disk regions 
can be heavily obscured by dust, and the sample completeness with 
respect to the parent stellar population light may also be affected. 

6 A stronger [Qui ] emission vs. Ha is f ound to be an excellent diagnostic 
tool for bright PNe 1 Ciardullo et al. 2004), while a reversed trend is likely to 
be expected for the fain test nebulae jMagrini et all200Cl) . See, for instance, 
Arnaboldi et al. 1 2003 ) and Corradi et al. 1 2005 ) for further diagnostic plots 
in the emission-line narrow-band color domain. 
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Figure 9. The expected bolometric correction for theoretical template 
galaxy models according to Buzzoni 1 2005). The models for different mor- 
phological type span the age range from 1 to 15 Gyr (the latter limit being 
marked by the solid dot on each curve). Bolometric correction refers to the 
V (upper panel) and B band (lower panel). A value of (Bol — V) = 
BCv = —0.85 mag can be taken as a representative correction for the 
whole galaxy types within a 10% uncertainty, as shown by the arrow on the 
upper plot. This also translates into BCb — —0.85 — (B — V) ga i for the 
B-band correction, as displayed by the dashed line in the lower panel. 



A comparison of the ljacobd fcOOfJ) updated SMC PNLF with 
the standard double-exponential function, as in eq. 1 1 81 . is proposed 
in Tableland Fig. [5] If the faint-end tail of the PNLF indeed oc- 
curs at 8 magnitudes down from the bright cut-off, then one can 
determine the fraction of PNe in the brightest 2.5 mag range, and 
define a parameter 02.5 as the luminosity-specific number of PNe 
with Al — Al* ^ 2.5. Figure[S|and Table|5|show that such bright 
nebulae represent a fraction between 10% and 20% of the total PN 
population. 

As we will see in the following sections, the definition of 02.5, 
which includes only the brighter PNe, makes the comparison with 
observations in external galaxies easier. For our discussion, we as- 
sume a simple relation such as a = 10 x 02.5, according to the 
standard PNLF (column 4 in Table |3j- This normalization is cor- 
rect if the double-exponential PNLF formula applies, and there are 
no PNe 8 mag fainter than the bright cut-off. 



3.1.1 Galaxy bolometric correction 

To compare observations and model predictions for the luminosity- 
specific PN number in external galaxies we need to convert the 
monochromatic galaxy luminosity to bolometric. One can use the 
Johnson V or B magnitudes for galaxy photometry, corrected for 
the distance modulus such as to match the absolute scale, and usu- 
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Table 4. The PN census in the Local Group galaxies'") 



Name 


Morph. 
Type 


Distance 
Mpe 


M B 


(S - V) 




Observed 
no. of PNe 


Completeness limit 
M lim -M* 


Standard 


ot 

SMC 


log a ±ff 


Reference 


M31 (all) 


Sb 


0.76 


-21.55 


0.68 ± 0.02 


+0.1 


~2700 




19000 ± 8000 




(i) 


M31 (bulge) 






-18.01 


0.95 ± 0.02< d ) 




94 


2.5 


940 ± 97 


497 ± 51 


-6.94g;l| 


(2) 


Milky Way 


Sbc 


0.01 


-20.80 


0.63 l -"'> 


-0.2 


~2000 




25000 ± 19000 


-6.408:11 


(3), (4) 


M33 


Scd 


0.80 


-18.74 


0.47 ± 0.02 


-0.5 


152 


2.6 


765 ± 85 


415 ± 46 


7 13 014 




LMC 


SBm 


0.050 


-17.93 


0.44 ± 0.03 


-0.52 


~1000 
~350 


~ 7.0 
~ 5.0 


1040 ± 60 


960 ± 50 


D - D '0.04 


(5) 
(6) 


SMC 


SBm 


0.060 


-16.24 


0.41 ± 0.03 


-0.84 


105 


6.0 




167 ± 19 


-fi R70-05 
D -°'0.05 


(7),(6) 


NGC205 


E5 


0.76 


-15.96 


0.82 ± 0.05 


-0.3 


35 


3.5 


134 ± 28 


87 ± 18 


-6 - 88 0\22 




M32 


El 


0.76 


-15.93 


0.88 ± 0.01 


-0.58 


30 


2.4 


186 ± 44 


97 ± 23 


' ' '0.31 




IC10 


Im 


0.66 


-15.57 


0.58 ± 0.05 


-0.7 


16 


2.0 


153 ± 46 


72 ± 22 


_ R cn0.20 
D - Dy 0.40 




NGC6822 


to 


0.50 


-15.22 


0.47 ± 0.15< c > 


-0.62 


17 


3.5 


74 ± 21 


48 ± 13 


-6.69g;lf 




NGC 185 


E3 


0.66 


-14.90 


0.73 ± 0.01 


-1.0 


5 


2.7 


45 ± 20 


25 ± 11 






NGC 147 


E5 


0.66 


-14.48 


0.78 ± 0.05 




9 


3.9 


29 ± 12 


20 ± 8 


-fi 91 019 




Sex A 


to 


0.86 


-13.02 


0.37 ± 0.08 


-1.32 


1 


1.9 


16 ± 16 


7 ±7 


-fi.38^ 30 




SexB 


to 


0.86 


-12.96 


0.51 ± 0.03 


-0.75 


5 


3.0 


37 ± 17 


22 ± 10 


-fi in - 21 




Leo A 


to 


0.69 


-11.36 


0.31 ± 0.08 


-1.51 


1 


3.0 


8 ± 8 


5 ±5 


-S.99^ 30 





References: (1) Merrett et al 2003 1; (2) Ciardullo et al. 1989); (3) Jacobv 1980; (4) Alloin et al. 1976 1; (5) Reid & Parker, 

I2002l l 



private communication: 



; f6) lJacob\ll200fll : (7) lJacobv & De Marcol 



( a ) Gala xy morphological types an d (B — V) colors are from the RC 3 catalog de Vaucoiileurs et a! 1991 ), distances and metallicities 

from Corradi & Magrini 2006). absolute B magnitudes, M s , from Karachentsev 2005) rescaled to the Corradi & Magrini 2006 ) distance modulus. Both (B — V) and M E 
quantities are reddening- and inclination-corrected according to the corresponding literature sources. 

< k > Adopted PN Oxygen abundance, [O / H] = log(0/H) - log(O/H) . For the Sun, 12 + log(O/-H) = 8.87 iGrevesse et alll996ll . 

( c ) PN population size within 8 mag from the PNLF bright-end tail, inferred according to two different extrapolation methods: 

Standard: empirical PNLF from Jacobv 1989), as in Tablel3Tcolumns 2 and 4), corrected for metallicity after Dopita et al. 1992) 

SMC: using the Jacobv 200fij) observed SMC PNLF (see columns 3 and 5 in Tablel3l. and correcting for metallicity as for the "Standard" case. 

Data from HyperLeda Lyon. 

( e ) The (B — V) estimate for the Milky Way is from the Robin et al. 2003 ) synthetic model, as quoted by Boissier & Prantzos 1999). 



ally complemented with a (rough) estimate of the color excess, 
E(B — V), to account for Galactic reddening. 

Given the (reddening-corrected) values of Mb and My, the 
required transformations to provide L ga i in eq. I16i are the follow- 
ing: 



IQ-0A (Afe-5.41) lQ-0 i (BCb+0.69) 
■^q-0.4 (My -4.79) 1Q-0.4 (BCV+0.07) 



(20) 



where BCb and BCv are the bolometric corrections to the B and 
V band, respectively; according to lB89t in our notation the Sun 
has M£ o1 = +4.72 mag, BC% = -0.07 mag, and BCq = 
—0.69 mag. A direct estimate of the bolometric correction from 
the galaxy observations is not a straightforward task, and the alter- 
native way is to rely on models. 

Figure[9]shows the trend of the B and V bolometric correction 
for the lBuzzonl ^2005) galaxy templates over the age range from 1 
to 15 Gyr. One sees that a simple and convenient solution, within a 
10% internal accuracy (i.e. ±0.1 mag) in the transformation, can be 
found for the V-band correction, assuming a fixed value BCv = 
(Bol — V) — —0.85 mag all over the relevant range of galaxy 
types and age. A suitable estimate for the B-correction is BCb = 
(Bol - V)- (B-V) = -0.85 — (B — V) mag. 



4 COMPARISON WITH OBSERVED PN POPULATIONS 
4.1 The Local Group PN census 

The LG galaxies have been extensively search ed for PNe since the 
discovery of five PNe in M31 by Baade ( 1955) . Recent surveys in- 
cluding all LG galaxies with log(Ls / Lq) > 6.7, provide a com- 
prehensive view of LG PN population which can be used as a first 
comparison for our population synthesis models. An updated list of 
the number of PNe known in the LG can be found in the work of 
ICorradi & Magrinil l200d) . While PNe have been found in 20 LG 
galaxies, the survey completeness and photometric accuracy are not 
good enough to allow a confident estimate of their global popula- 
tion in all of them. 

In Table |4] we therefore summarize the updated information 
for the PN census only for those LG galaxies whose survey com- 
pleteness limit (and the number of PNe within it) is known. Many 
of the d ata come from the so-called Local Group Census project 
(see e.g. lCorradi et alj 2005. and references therein); other sources 
are indicated in the table. 

From these data we have estimated the total PN population 
size of these galaxies (-/V to t in Table |4j in two ways, to account 
for the uncertainty in our knowledge of the shape of the PNLF at 
faint magnitudes. First, for each galaxy the magnitude difference 
between the completeness limit of the survey and the expected ap- 
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Figure 10. A comprehensive overview of the luminosity-specific PN number in LG galaxies (star markers) and external ellipticals from Tablel5l(solid triangles) 
and Table lol (open triangles). PN d ata for local gala xies are from Table l4l and are based on the "Local Group Census Project" of lCorradi et alJ 120051) . Also 
superposed on the plot, there are the Bu zzoniH2005l) template galaxy models, as summarized in Tablel2l Galaxy evolution is tracked by models along the whole 
E-Sa-Sb-Sc-Sd-Im Hubble morphological sequence from 5 to 15 Gyr, with the latter limit marked by the big solid dots. Two model sequences are reported on 
the plot assuming an IFM R as from the standard case of a Reimers mass loss parameter rj = 0.3 (lower sequence), and from the empirical relation of lWOfJ 
(upper sequence). For the WOO models, the relevant data of Tablel2lhave been corrected by A.(B — V) = —0.02 mag and log a = log aw + 0.04, according 
to the arguments of footnote J5J- An indicative estimate of the mean representative PN lifetime (in years) is sketched on the right scale, according to eq. 1221 . 



parent magnitude of the PNLF cut-off, m*, was computed. The 
value of m* is determined assuming M* = —4.47, the metallicity 
correction in eq. 1 1 91 with the value of [O/H] as report ed in column 
6 of the table, and the distances as in ICorra di & Maerini 1 2006) 
(listed again in column 3 here). The observed number of PNe within 
the completeness limit has then been extrapolated 8 mag down the 
PNLF cut-off using the empirical formula in eq. 1 1 8i . as reported 
in Table|5J resulting figures are listed in column 9 of Table|4| The 
total PN populatio n size was also estimated rescaling the observed 
PNLF of the SMC Ijacobv & De Marcol20oi see column 5 of Ta- 
ble[3j, which is complete 6 mag down the cut-off, and assuming a 
~ 50 % incompleteness for the next magnitude bin, consistent with 
recent deeper observations I Jacobv 2006). The corresponding total 
PN population estimated for the LG galaxies is indicated in column 
10 of Table|4| 

For the S MC, the adopted PN population is the number es- 
timated by I Jacobv & De Marcel |2002), plus some 20% to includ e 
the faintest PNe as suggested by recent observations I Jacobv 2006). 
The LMC population size is determined using approximate figures 
of the observed numb er of PNe an d the depth of the discovery sur- 
veys, according to Jacobv 1 2006 ). For M 31, data for a relatively 
dust-free region in the bulge iCiardullo et all l989) were adopted; 
the local number of PNe was then rescaled to the whole galaxy lu- 
minosity. The latter value represents an estimate of the global pop- 
ulation, and the derived value of a reported in Table [4] is in fact 
representative of the bulge evolutionary environment. Finally, the 
total PN population for the Milky Way in Table|4]is a conservative 
estimate from the data bv lJacobvl ll98Ch and lAlloin et alJ Jl976l) . 

The difference between the total PN population based either 
on the "standard" or the "SMC" PNLF is of a factor two or smaller. 



The derived value of a is reported in column 1 1 of Table [4] as- 
sumin g for LG galaxies the B — V colors from the RC3 cata- 
log ide Vaucouleurs et al. 1991) and the absolute Mb magnitudes 
from Karachentsev 1 2005), converted to bolometric according to 
eq. I20i . 7 The quoted error bars in the Table are conservative esti- 
mates of the uncertainty of iVtot , which account for the difference 
between the "standard" and "SMC" PNLF extrapolated values. 

The LG galaxy sample provides the natural benchmark to test 
our theoretical models over a range of evolutiona ry environ ments 
and star formation histories. Our predictions for the B uzzon il 120051) 
template galaxy models of Tableware compared in Fig.^|with the 
LG data of Table|4]for a relatively old age range, from 5 to 15 Gyr, 
and mass loss scenarios (rj = 0.3 and WOO). 

The remarkable feature of the log a plot, for the observed 
galaxies representative of the whole Hubble morphological se- 
quence, is the fairly constant PN rate per unit galaxy luminosity. 
Data support an average rate between 1 and 6 PNe per 10 7 Lq. 
Such a value is related to a narrow range of PAGB stellar core 
mass, according to the calibration of Fig. [5] with M core less than 
0.60-0.65 M@. 

There are at least three important consequences of this sharp 
mass distribution: i) the mass-loss scenario supported by the obser- 
vations better agrees with the WOO IFMR, which implies a stronger 

7 To be consistent with theJCorrad^&MagrinJ l2006t) survey and com- 
pleteness correction, the Karachent se\n2005ir*absolute M b magnitudes 
have been slightly rescaled to the Corradi & Maarini 1 2006) adopted galaxy 
distances, as reported in Table |4] Note, in addition, that both Mb and 
(B — V) are reddening- and inclination-corrected values according to the 
original data sources. 



12 A. Buzzoni, M. Arnaboldi & R.L.M. Corradi 

Table 5. Recent additions to PN census in early-type galaxies 



Name 



Morph. 
TypeM 



Observed Comp. limit 



N t 



(b) 



log a 



Ref. 



no. of PNe 



M-M* 



NGC 1316 


SO 


43 


1.0 


1720 ± 262 


-7.50 ±0.07 


Arnaboldi et al. (1998) 


NGC 1344 


E5 


197 


1.0 


2300 ± 47 


-6.75 ±0.02 


Teodorescu et al. (2006) 


NGC 1399 


El 


37 


1.0 


1480 ± 243 


-7.30 ±0.07 


Arnaboldi et al. (1994) 


NGC 3115 


SO 


61 


1.0 


2440 ± 312 


-6.59 ±0.05 


Ciardullo et al. (2002a) 


NGC 3379 


El 


109 


2.0 


1535 ± 104 


-6.77 ±0.03 


Romanowskv et al. (2003) 


NGC 4697 


E6 


535 


2.5 


3500 ±231 


-6.82 ±0.03 


Mendez et al. (2001) 


NGC 5128 


SO 


431 


2.5 


4300 ± 207 


-6.30 ±0.02 


Hui et al. (1993) 



( a ) From the RC3 catalog <de Vaucouleurs et alll99lh . 

(fc) PN population size within 8 mag from the PNLF bright-end cut-off, inferred according to standard 
PNLF, as in Tablel3*l(columns 2 and 4). 



mass loss for intermediate and high-mass stars compared to the 
standard scenario (77 ~ 0.3-0.5) for Pop II stars as in Galactic 
globular clusters; ii) according to Sec. 12. 21 for a consistent fraction 
of the PN population in LG galaxies the inferred lifetime is con- 
strained by the dynamical timescale of nebula evaporation rather 
than the stellar core mass evolution; Hi) the latter evidence also 
supports a small dependence of a with time and distance (see right 
panel of Fig. pointing to a relatively "universal" shape for the 
PNLF (but see, however, Sec. l5.2l for an important warning in this 
regard). 

4.2 PN surveys in external galaxies 

As discussed in Sec. 13. II a complete survey of the PN population 
in late-type galaxies is often plagued by spurious detections caused 
by Hll regions, SN remnants and, for galaxies at distances larger 
than 10 Mpc, by Lya background galaxies at z ~ 3.13, (see e.g. 
ICiardullo et all2002dlArnaboldl2004l) . 

Current surveys only investigate the brightest magnitude bin 
of the PNLF, to measure the bright cut-off magnitu de M* and 
lead therefrom to galaxy distance (see, e.g. lCiardullo et ai] [2002a. 
for a recent survey of six SO and active, mostly Seyfert 2, spi- 
rals). The best samples o f PNe in late-type systems still remain the 
ones for M31 and M 33 jMerrett et alJl2003t iMagrini et alj|200fi 
ICiardullo et all2004) . 



4.2.1 PN samples in early-type galaxies 

The lack of active star formation and a negligible fraction of resid- 
ual gas in early-type galaxies make the [Om] PN detection rel- 
atively straightforward across the whole galaxy body, excluding 
the innermost regions, where the galaxy spectral continuum is too 
bright. 

However, only low-mass ellipticals and dwarf spheroidals can 
be found among the LG galaxy population and, with the exception 
of NGC 5128 (Cen A) at 3.5 Mpc, the nearest normal and/or gi- 
ant ellipticals are at 10 Mpc or further. At these distances, the 
typical [Om] magnitudes are in the range 26-27 mag or fainter, 
and spectroscopic-confirmed PN samples are usually limited to 
M < M* + 1.0, with only few good cases observed in slitless 



spectroscopy. Hopefully, larger PN samples in early-type galaxies 
will soon be available with the dedicated PN.S spectrograph oper- 
ating at the William Herschel Telescope at La Palma I Dougl as"et alJ 



A comprehensive collection of PN data for dwarfs, giant el- 
liptical and SO galaxies in the LG , Leo group and the Virgo cluster 
can be found in lHuietalJiT99l . comparing original observations 
of NGC 5128 with the value of 02.5 for a sample of 13 early-type 
galaxies plus the bulge of M 31. An updated census for some of 
these objects, plus additional results from deeper surveys for a few 
more galaxies are listed in Table|5] In column 4 we report the total 
number of detected PNe (at all magnitudes) together with the com- 
pleteness limit of the survey, (M — A/*), and the inferred number 
of the global PN population on the sampled galaxy region (with its 
Poissonian error estimate), by assuming a standard PNLF down to 
M = M* + 8.0 mag. The value of a is then computed by normal- 
izing to the sampled galaxy luminosity, according to the original 
data sources in the literature. 



For the data of Table [5] and the lHui et alJ (l993) sample, we 
collected, in Table [6] supplementary dynamical and photometric 
information from iBurstein et alJ <1988l) . the RC3, NED and Hy- 
perLeda on-line databases. Two LG galaxies were also added to 
the sample, including the dwarf satellite ellipticals of M 31 (i.e. 
NGC 205 and M 32), and the bulge data for M 3 1 itself. This whole 
sample of early-type gala xies is displayed in Fig. 1101 matching the 
LG data and the Buzzoni 1 2005) template galaxy models. 



As pointed out bv lHui et alJ 11993T) . in the plot of Fig. llOl one 
can notice a sharp decrease of the luminosity-specific PN number in 
early-type galaxies compared with the theoretical predictions of our 
E-galaxy model and the empirical estimate of a for the Milk y Wa; 
and o ther nearby sp irals. In addition, following Peimbert 1 199' 
and H ui et alJll993l) . the figure also reports a clear trend of a with 
galaxy color, with a poorer PN population (per unit bolometric lu- 
minosity) in redder ellipticals. So far, this trend has not received a 
satisfactory explanation. 
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Table 6. The updated PN sample for local and distant early-type galaxies^' 



NGC 


logo- 


(B - V)o 


Mg 2 


(1550 - V) 


logaW 


Ref. 


Notes 




[kins" 1 ] 


[mag] 


[mag] 


[mag] 








205 


1.61 


0.82 


0.071 


1.19 


-6.88 


(1) 


star forming 


221 


1.90 


0.88 


0.198 


4.50 


-6.77 


(1) 


M32 


224 


2.27 


0.95 


0.324 


3.51 


-6.94 


(1) 


M31 (bulge only) 


1316( c > 


2.38 


0.87 


0.265 


5.0 ( d > 


-7.50 


(2) 


For A - merger 


1344 


2.22 


0.85 


0.267 




-6.75 


(2) 




1399 


2.52 


0.95 


0.357 


2.05 


-7.30 


(2) 




303l( c ) 


2.23 


0.82 


0.295 




-6.80 


(1) 




3115 


2.45 


0.94 


0.309 


3.43 


-6.59 


(2) 




3377' c ) 


2.16 


0.84 


0.273 




-6.43 


(1) 




3379 


2.33 


0.94 


0.329 


3.86 


-6.77 


(2) 




3384< c ) 


2.20 


0.91 


0.296 


3.9 < d > 


-6.42 


(1) 




4374 


2.48 


0.94 


0.323 


3.55 


-6.77 


(1) 




4382 


2.24 


0.88 


0.242 


4.22 


-6.70 


(1) 




4406 


2.42 


0.90 


0.330 


3.72 


-6.89 


(1) 




4472 


2.49 


0.95 


0.331 


3.42 


-7.16 


(1) 




4486 


2.60 


0.93 


0.303 


2.04 


-7.10 


(1) 




4594(c) 


2.41 


0.84 


0.340 




-6.66 


(1) 




4649 


2.56 


0.95 


0.360 


2.24 


-7.22 


(1) 




4697 


2.25 


0.89 


0.320 


3.41 


-6.82 


(2) 




5128M 


2.14 


0.88 




> 5 ?< e ) 


-6.30 


(1) 


Cen A - merger 



References for a estimates: (1) Hui et al, 1 1993); (2) this paper 



Stellar velocity dispersion log a, Mg2 index and (1550 — V) color from Burstein et al. 1 1988) unless 
other wise stated ; redde ning-corrected (B — V) from the RC3 catalog ideVaucoTleur^taiTl99ll) ; 
W For the lHuietal]<1993l) data, log a = log o 2 .5 + 1-00. 
' c ) log <r and Mg2 from HyperLeda Lyon; 

( ' (1550 — V) color estimated from the ultraviolet galaxy catalog o f lRifatto etajll 19951) and RC3/HyperLeda 
dereddened total V magnitud e. Reddening correction is accordingTo Tseatorni979l) . assuming 
a Galaxy extinction map from Burstein & Heiles 11984) : 

' e ' Conservative estimate on the basis of the upper limit to the 1540 A galaxy flux, as reported by the NED 
database. 



5 THE PN POPULATION IN ELLIPTICALS 

From the definition of a in eq. a lower value of the luminosity- 
specific PN number corresponds to a shorter PN lifetime. If 
30 000 yr is a reasonable value for tpn in the most PN-rich galax- 
ies, then 



For a low value of a, like for instance in NGC 13 16 or NGC 1399, 
one leads therefore to a lifetime of 3000-5000 yr for the galaxy PN 
population (see the right-scale calibration in Fig. I10i . 8 The value 
of tpn linked to the observed value of a in eq. 12 1 1 must, how- 
ever, be regarded as an average on the PAGB core mass distribu- 
tion, which cannot be determined observationally in systems other 
than the Milky Way (e.g. lZhang & Kwokll993l) . 

As far as we restrain to standard stellar evolution theory, such 
a shorter (or even a vanishing) PN lifetime can be obtained in three 
different (and to some extent mutually exclusive) ways, either 

(i) by increasing the stellar core mass so that nuclear evolution 
speeds up (see footnote^, 



8 Note that our definition of the PN lifetime relates to the nebula visibility 
phase, and it is different from the kinematic age of the system, as derived 
from the ratio between a bsolute size of the m ain nebular shell and its ex- 
pansion velocity (see, e.g. Villaver et al. 2002, for a discussion). 



(ii) by delaying the Hot-PAGB phase of the stellar core (for in- 
stance, by increasing Tit in eq.[5) such as to let the nebula evaporate 
before it can be excited, or 

(iii) by fully inhibiting the AGB phase so that the PN event can- 
not take place, at least in a fraction of the galaxy stellar population. 

Each of these different scenarios leaves a different "signature" in 
the overall shape of the PNLF and its bright cut-off luminosity. 
Definitely, case (i) above is the less favoured one in our analysis, 
as it would predict bluer ellipticals to have a lower value of a, con- 
trary to what we observe. Likely, case ( ii) and ( iii) better fit with an 
overall evolutionary scenario dominated by low-mass (old?) PAGB 
stars, and they could be at work at the same time (though to a differ- 
ent relative degree) among the PN population of early-type galax- 
ies, as the recent HST observations of M 32 iBrown et alJl200oT) 
seem to support. According to the arguments of Sec. 12.3.31 the core 
mass range 0.52 < M cola < 0.55 Mq may be the preferred one 
to explain the case (ii) scenario, while a more intense core mass 
depletion, leading to Af corc < 0.52 Mq, would cause a fraction of 
Post-HB stars to override the PN event, like in case (iii). 

5.1 PN core mass and M* invariance 

Within some still larg e theoretical uncertainties, stellar evolution 
models iMendez & SoffneJl997i:lMarigo et all2004l) a gree on the 
fact that the brightest PNe are not always related to the most mas- 
sive cores. Nonetheless, stars of Mi > 2Mq, ending up their 
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Figure 11. The observed distribution of the elliptical galaxy sample of Ta- 
ble |6| (plus M 32 and the bulge of M 31, as labeled on the plot) versus 
Lick spectrophotometric index Mg2 . Note the relative lack of PNe (per unit 
galaxy luminosity) in more metal rich ellipticals. The relevant case of the 
merger galaxy NGC 1316 is singled out, while the two active star forming 
ellipticals NGC 205 and NGC 5102 are out of range with Mg 2 < 0.1 and 
not shown. See text for a discussion. 

PAGB evolution with M corc > 0.7 M (cf. Fig.Et, are required to 
generate the M* nebulae (see, e.g. Fig. 10 in M arigo et alj2 004). 

If our models for PN evolution are correct, then one predicts 
low-mass nebulae (M cole < 0.65 Mq) to dominate the PNLF of 
galaxies of different morphological type. In fact, such a claim is 
even stronger for ellipticals, for which very low core-mass values 
must be invoked on average for their old metal-rich stellar pop- 
ulations. Therefore, it becomes harder in these galaxies to justify 
the presence of relatively high-mass nebulae reaching the M* lu- 
minosity, as the empirical invariance of the PNLF bright cut-off 
magnitude may imply. 

A prag matic approach to the problem has recently been 
pursued by ICiardullo et alj 120051) . who argued on the possi- 
ble presence of a blue-straggler (BS) stellar population in ellip- 
ticals, via coalescence of close binary systems. BSs are com- 
monly observed in Galactic ope n clusters of all ages (e.g. lKinmanl 
1965t lEggen & Sandagd [l969 ) and in some globulars, too (e.g. 
BTionarmoeHu^l994l) . and these objects may in principle reach up 
to twice the TO mass, that is ~ 2 Mq, even in old stellar systems. 
While ensuring a convenient fraction of M* nebulae and account 
for the observed M* invariance, this scenario still leaves, however, 
a few open question, which we address below. 

5.2 Diagnostic planes 

To evaluate the impact of the different AGB and PAGB evolution- 
ary pictures in elliptical galaxies, it may be useful to investigate 
the correlations of a with other observed quantities for the galaxy 
sample of Table|6| as displayed in Fig. ll llll2l and|13| 

The plot of the Lick Mg 2 index jFaber et all 19851) in Fig.fTTI 
is of special interest, in this regard, because it is the least affected 
by reddening. Here, the a correlation with galaxy color is more 
cleanly replicated, suggesting that metallicity, rather than age, is 
the relevant parameter that affects the observed PN rate per unit 
galaxy luminosity. With the exception of NGC 1316 (Fornax A), 
that we shall discuss in some details, the plot shows that a lower 



Figure 12. Same as Fig. II II but for the galaxy velocity dispersion a in 
km s — . It is evident a lower value of a in high-a (roughly more massive) 
galaxies. See text for further details. 

number of PNe per unit galaxy luminosity is produced in metal- 
rich ellipticals. 

As chemical and dynamical prop erties are tig htly correlated in 
early-type galaxies iFaber & Jacksor]ll976t iTerlevich et aljfl98ll 
iBurstein et alJll984T) . one may expect some correlation of a with 
the galaxy internal velocity dispersion too (that is, roughly, with 
the galaxy total mass). This is shown in Fig. 1121 Among others, 
the good fit of the M 3 1 bulge to the overall correlation for early- 
type systems in Fig. ll ll and !12l could be considered as an additional 
piece of evidence of the similarity between the stellar populations 
of spiral bulges and ellipticals l Jablonka et al. 199a). 

The general picture of PN evolution, sketched in Sec. 12.31 
especially reflects in the distinctive properties of the early-type 
galaxy population in the ultraviolet spectral range. In particular, 
the two relevant Post-HB evolutionary paths, that lead stars ei- 
ther to a full AGB completion or straight to the high-temperature 
white-dwarf cooling sequence have different impact on the galaxy 
spectral energy distribution (SED) at short wavelength compared 
to the optical luminosity. To explore this important feature, in 
Fig. ll3l we plot the log a values vs. the UV color ( 1550 — V) = 
-2.5 log[/(1550 A)//(V)], as first defined by IBurstein etafl 
( 1988), where the galaxy SED is measured at 1550 A and in the 
Johnson V band. The (1550 — V) color, from IUE observations 
and co rrected for Galaxy extinction, was provided by Burstein et alJ 
( 1988) for a fraction of galaxies in Table|6| 9 

Figure ^| shows a tight correlation between the value of a 
and the ultraviolet emission, with massive UV-enhanced ellipticals, 
like NGC 4649 and NGC 4486 in Virgo, or NGC 1399 in Fornax 
that are also much poorer in PNe. To a finer analysis of the plot, one 
can even notice an apparent gap (marginally evident also in Fig. 1121 
between these three giant ellipticals and the bulk of more "normal" 
galaxies (including the bulge of M 3 1). 

One may speculate that the bulk of the PN population in "nor- 

9 The (1550 - V) color for NGC 13 16 and NGC 3384 h as also been 
added to Tabled] as obtained from the iRifatto et al] 1 1995) UV catalog. 
These entries, however, are reported with a larger error, as the galaxy flux 
at 1550 A derives from a crude extrapolation of the reported magnitudes at 
1650 A and 2500 A. 
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Figure 13. The luminosity-specific PN numbe r versus u ltraviolet color 
(1550 — V), as originally defined by Burstein et al. 1 1988), for the ellipti- 
cal galaxy sample of Tablel6*l(plus the Andromeda satellites and the bulge 
of M 31). Some relevant cases, like NGC 205 and NGC 5102 (star form- 
ing), NGC 1316 and NGC 5128 (merger ellipticals) are singled out on the 
plot. Note the tight relationship between "quiescent" ellipticals and a, with 
UV-bright galaxies to be also PN-poor. See text for a full discussion of this 
important effect. 



mal" ellipticals evolves according to the case (ii) scenario; a sub- 
stantial fraction of stars, in these galaxies, would therefore com- 
plete its AGB evolution leading to low-mass PNe, with M corc ^ 
0.55 Mq . The intervening increase of the AGB =>■ Hot-PAGB tran- 
sition time in this mass range, would eventually shorten the nebula 
lifetime and reduce tpn from ~ 30 000 yr to ~ 10 000 yr or less, 
as for the stellar population in the M 31 bulge (see Fig. llOi . If the 
mean PAGB core mass is further decreased in more massive galax- 
ies, perhaps as a consequence of a stronger mass loss, then one may 
expect an increasing fraction of HB stars to feed the AGB-manque 
evolutionary channel and fail to produce PNe, as in case (Hi) sce- 
nario. This would imply an even lower value of a and a strongly 
enhanced galaxy ultraviolet emission, as observed in NGC 4649, 
NGC 4486 and NGC 1399 in Fig.fBl 

In this framework, the envisaged role of the BS stellar com- 
ponent could be assessed on the basis of the Kin & D eng 1 2005) 
analysis of a sample of old (t 5 Gyr) Galactic open clusters. 
A sizable presence of BSs is found to severely affect cluster col- 
ors, by shifting the integrated B — V over ~ 0.2 mag to the blue 
(see Fig. ll4i . When applied to ellipticals, this argument could place 
a quite tight constraint to the overall BS luminosity contribution, 
and the size of the induced PN progeny, in order to avoid unreal- 
istically bluer galaxy colors. For example, if we take the case of 
cluster M 67 (the Praesepe) as a refere nce, BSs are foun d to sup- 
ply ~ 50% of the total B luminosity <Deng etalJll999l) . causing 
a A(B — V) = — 0.15 mag shift to the integrated cluster color 
IXin & Deng||2005l) . Converting these values to bolometric lumi- 
nosities (e.g. by relying on our discussion in Sec. 13.1.11 . this im- 
plies that BS contribution must be well below ~ 10% of the total 
galaxy luminosity if we wish to avoid any measurable effect on the 
integrated optical colors of ellipticals. 

In addition, if M* PNe in ellipticals really stem from BS evo- 
lution, then one consequence is that the value of a for these galax- 
ies, as extrapolated from the observation of the very brightest nebu- 
lae alone, could sensibly overestimate the population of "standard" 



Figure 14. The blue shift of the integrated B — V color of old open clusters 
in the Galaxy caused by the BS stellar population, from Xin & Dena 1 2005 ). 
The BS component (N^s) is normalised in terms of its ratio to the number 
of MS stars down to 2 mag below the TO luminosity (N2). Solid dots are 
for the oldest (t 5 Gyr) clusters, while star markers include clusters with 
1 t < 5 Gyr. Symbol size is proportional to cluster statistical richness. 



(i.e. single-star) PNe of lower core mass. According to eq. 12 H . this 
would imply an even shorter mean lifetime for single-star nebulae, 
thus enforcing the importance of the case (ii) and (Hi) evolution- 
ary channels. If two PN sub-populations of bright (coalesced) PNe 
and fainter nebulae generated from standard evolution of low-mass 
stars really coexist in early-type galaxies, then one might wonder 
whether they also display any distinctive dynamical signature, as 
the recent case of NGC 4697 l Sambhu s et afll2006l) seem to sug- 
gest. 

To complete the discussion of Fig. 1131 we now focus on two 
major outliers in the plot, namely NGC 205 and NGC 1316. The 
first case can easily be assessed as this M 31 dwar f satellite is 
know n to undergo active star formation (e.g. lBica et all 1 99(1 iLed 
1 19961) and its strong ultraviolet emission (as well as its exceedingly 
low Mg2 index, see Table [6) is in fact the result of young (a few 
10 8 yr) MS stars. More special attention, instead, must be paid to 
the case of NGC 1316 (Fornax A). 



5.3 Galaxy environment and PN evolution: the case of 
NGC 1316 and Cen A 

Together with NGC 5128 (Cen A), NGC 1316 is known as one of 
the best established exampl es of "mergers" a mong early-type sys- 
tems <SchweizeJl983l : lGoudfrooii et all200ll) . The study of these 
two galaxies may therefore lead to a preliminary assessment of the 
influence of "active" galaxy environments on the PN population. 

The location of both objects in the log a plot of Fig. ll2| places 
them at the extreme (both highest and lowest) values of a; further- 
more, one must also report the peculiar location of NGC 1316 in 
the Mg2 plot (Fig. II II . which shows a severe PN deficiency when 
comparing for instance with NGC 4382, of similar mass and metal- 
licity. 10 Regarding the ultraviolet properties of these "merger tem- 
plates", the upper limit to the 1540 A flux for NGC 5128, as re- 



10 Unfortunately, the lack of Mg2 data for NGC 5128 prevents a similar 
comparison for this galaxy. 
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Figure 15. A representation of the envisaged PN evolution versus core mass of PAGB stars. The effect of different parameters, like metallicity, mass loss and 
age is outlined. In particular three evolutionary regimes are singled out, with PN visibility lifetime rp^ (and correspondingly a) constrained respectively by 
the nuclear timescale (tpagbX shell dynamics (r dyn ), and transition time (r t t). PN visibility drastically reduces for M core < 0.55 Mq until reaching a 
critical limit for PN formation about M core ~ 0.52 Mq. See text for full discussion. 



ported by the NED database, suggests a presumably very "red" UV 
color for this galaxy, such as (1550 — V) S> 5, in line with the 
observations of NGC 1316. 

Definitely, in spite of similar observed spectrophotometric 
properties (i.e log a, (B — V) and, likely, the (1550 — V) color, 
cf. Table©, it remains difficult to understand the role of galaxy 
merging mechanisms that led, for Cen A and For A, to such a dra- 
matically different behaviour in terms of PN population. 



6 DISCUSSION AND CONCLUSIONS 

Wide-field CCD detectors, and accurate selection criteria based on 
Ha and [Olll] narrow-band photometry, allow a systematic survey 
of the PN population both in the Milky Way and in external galax- 
ies. Deep observations now explore the PN luminosity function sev- 
eral magnitudes below the bright cut-off, as recently achieved in 
the SMC (Jacobv 2006). The corresponding developments in the 
theoretical modeling help in clarifying the main physical mecha- 
nisms that constrain the nebular properties, especially those related 
to shell dynamics fe.g. lVillaver et al|2002t IPerinotto et aljl2 004a) 
and chemical composition (Liu et al. 2004; Perinotto et al. 2004b), 
although a clear understanding of the the PNLF at bright and low 
magnitudes is still to come. 

One open question deals with the empirical evidence for a 
nearly constant absolute magnitude of the PNLF bright cut-off 
(M*), which makes bright PNe effective standard ca ndles for the 
intermediat e cosmic d istance scale (< 100 Mpc; Jacobvl ll989t 
ICiardullo et ai] E002a ) . This feature is not explained by the theo- 
retical models, which predict Al* to be a function of age, metallic- 
ity and other distinctive properties of the parent stellar population 
I Marigo et al. 2004). In particular, this problem is stronger for ellip- 
tical galaxies, where a lower TO mass for their stellar populations 
may hardly reconcile with the required pres ence of ~ 2 Mp, stellar 
proge nitors for the M * nebulae l Ciardullo et al. 2005 ; Marig oet alj 



The PNLF faint-end tail is also subject to a substantial uncer- 
tainty, with observations showing a decrease in th e number of PN 
at about 6 mag fainter than M * IJacobv & De Marcoll2002l) . and 
theory which predicts a distribution reaching magnitudes as faint as 
8 mag below the PNLF bright cut-off. We have discussed in Sec. 14. II 
that our limited knowledge of the PNLF faint end reflects in a factor 
of two uncertainty on the estimated total number of nebulae for a 
given stellar system. The errors become larger for the most distant 
galaxies, external to the LG. In most cases, the PNLF is only sam- 
pled in its brightest magnitude bin, and the observed PNe number 
must be multiplied by a factor of 20-40 (see Table|3j to provide an 
estimate for the whole PN population. 

Despite these large uncertainties, the study of the luminosity- 
specific PN number (the so-called "a" parameter in our discus- 
sion) in external galaxies has allowed us to tackle the properties 
of the PN population in different environments, from dense bulge- 
dominated galaxies to very low density stellar popula tions such as 
those of the intracluster dif fuse stellar component lAguerri et alJ 
l2005t Feldmeier et al. 2004). From the theory presented in Sec.[2| 
and [5] there is a close liason between the value of a and the PN 
visibility lifetime: as the specific evolutionary flux, B, in eq. |3} is 
nearly constant, then 

r PN « 30 000 (a/5.4 10~ 7 ) yr. (22) 

Depending on the stellar core mass at the beginning of the AGB 
phase, the value of tpm must be compared with three relevant 
timescales that determine the PN evolution. As sketched in Fig. 1151 
they are the PAGB core lifetime (tpagb), the dynamical timescale 
for the nebula evaporation (rdyn), and the transition time for the 
stellar core to leave the AGB and reach the high-temperature 
regime required to trigger the nebula [OIII] emission (rtt). Along 
the SSP evolution, mass loss eventually settles the absolute clock 
that links the tpn evolution with the SSP age (see Sec. l2,4t . 

As far as core evolution provides the leading timescale for the 
PN visibility, standard mass loss theory a la Re imersl ll975h pre- 
dicts that t pn (and correspondingly a) should increase with SSP 
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age. 11 More generally, one would expect this to be the case of late- 
type galaxies, where star formation is on-going and a consistent 
fraction of high-mass stars is younger than a few Gyr (cf. Fig. [(J. 
As the fraction of young stars is higher in irregulars than in ellipti- 
cal galaxies, and the integrated B — V color becomes bluer along 
the E — » Im Hubble morphological sequence, then one may con- 
clude that lower values of a are expected in bluer galaxies. 

However, observationally this does not occur for the LG 
galaxy population, as shown in Fig. 1101 In fact, while the sug- 
gested theoretical range for a in the rj = 0.3 models matches the 
observations on average, the latter points to a constant behaviour 
with galaxy morphological type, with a typical rate between 1 and 
6 PNe per 10 7 L©. From these observations, the inferred PN life- 
time in LG spirals and irregulars should exceed 10000 yr, a value 
that requires the presence of a substantial fraction of stars with 
Mcore ;$ 0.65 Mq (cf. Fig. |5J even in those galaxies with the 
most active star formation activity. 

The relatively low final mass and the tight mass range required 
by the "dying" stars in external galaxies find independent confirma- 
tion in the Milky Way, as also indicated by the observational evi- 
dence of the IFMR for local open clusters by WOO. We have shown 
in Sec. 14. II that population synthesis models using this empirical 
IFMR produce a better fit (both in absolute value and relative trend 
with galaxy morphological type) for the correct value of PN density 
per unit galaxy luminosity in the LG galaxies, along the whole Hub- 
ble morphological sequence (see Fig. llOI . This result indicates that 
the dynamical evolution plays a central role in setting the overall 
PN observed properties as for low core-mass stars, the dynamical 
rather than nuclear timescale is the real driving parameter to con- 
strain PN visibility (see Fig. ll5> . As a consequence, one may state 
that, rather than probing any real mass distribution of PAGE stars, 
the PNLF basically tracks the time evolution of the expanding shell 
aroun d stellar nuclei of nearly fixed mass iHenize & Westerlundl 
Il963h . 

Within this picture, one problem is related with the PN evo- 
lution in (giant) elliptical galaxies. Here, the empirical evidence 
for a constant PNLF bright cut-off magnitude would require a siz- 
able component of relatively high-mass stars, well above the TO 
mass of ~ 1 Mq ex pected for these old stellar systems. Fol- 
lowing lciardullo et alJ l2005h . BS stars originating from coalesced 
close binary systems according to a classical evolutionary scheme 
lMcCreJl 964). may possibly overcome the dilemma. However, if 
ellipticals do indeed host a BS population in a fraction similar to 
what we detect in old Galactic open clusters, then galaxies should 
appear roughly 0.2 mag bluer than observed, in B — V (see Fig. 1141 . 
As a consequence, from our arguments in Sec. 15.21 one is led to 
conclude that, in any case, BSs cannot provide much more than a 
few percent of galaxy bolometric luminosity and PNe coming from 
the BS progeny must be confined to the very brightest bin of the 
PNLF, thus representing a marginal fraction of the global PN pop- 
ulation. 

The relative lack of massive stars in early-type galaxies, and 
a high (super-solar?) metallicity possibly easing a stronger mass 
loss in these galaxies, might actually lead to a larger component of 
blue HB stars, eventually evolving into low-mass PNe. In Sec. 12. 31 



11 Fo r a characteristic value of the mass loss parameter r\, we have from 
IR83 that the final mass of PN nuclei roughly scales as Mpagb oc 
^—0.35 ^—0.29 ^ reca jj m g that, to a first approximation, t oc M^Tq' 5 

lBuzzoni2003) . According to footnote^ this eventually leads to tpagb 
„2.2 t 1.8 



we have seen that for the low-mass cores (M CO rc <■ 0.55 Mq), 
the increase in the transition time r t t reduces effectively the visi- 
bility timescale of the nebula (see, again Fig. ll5i . thus decreasing 
the value of a. This trend is expected to depend on galaxy color 
(or metallicity, as shown by the Mg2 distribution of Fig. fTTl with a 
lower PN density per unit galaxy luminosity for the redder ellipti- 
cals. 

Furthermore, for the stellar component with M corc < 
0.52 Mq, both HB and Post-HB evolution occur at high tem- 
perature (T e ff ^> 10 4 K) thus skipping the PN phase e ntirely 
icastellani & Tornambeli99lllDorman et alll993tlB16ckeill995h . 
For these stars, the AGB manque evolution can effectively transfer 
some fraction of the PAGB energy budget to the hot HB tail and 
affect the integrated galaxy SED. As a consequence, one may ex- 
pect that the most UV-enhanced ellipticals display also the lowest 
values of a. We have shown in Sec.[5]that such a tight correlation is 
observed for elliptical galaxies in the Virgo cluster and other groups 
(see Fig. ll3i . Recent HST observations of the resolved stell ar pop- 
ulation of M 32 and the bulge of M 31 <Brown et alJl200Clll998l) 
further support this proposed scenario. 

The presence in our sample of two merger galaxies (For- 
nax A and Cen A) may provide the opportunity to study the effect 
of galaxy interactions on the PN population. Based on the avail- 
able data however, no firm conclusions can be drawn, although 
one has to remark that these two galaxies stand out as those with 
the most extreme (both highest and lowest) values of a in our 
sample. Perhaps such a huge variation in the galaxy PN popu- 
lation may be related to the disruptive effect of the i ntergalactic 
ram pressure, as recent l y co nsidered bv lVillaver et alJ l2003h and 
IVillaver & Stanghellinil 120051) . 
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